Experimental Study of the BEC-BCS Crossover Region in Lithium 6 
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We report Bose-Einstein condensation of weakly bound 6 Li2 molecules in a crossed optical trap 
near a Feshbach resonance. We measure a molecule-molecule scattering length of 170lgg° nm at 
770 G, in good agreement with theory. We study the 2D expansion of the cloud and show deviation 
from hydrodynamic behavior in the BEC-BCS crossover region. 

PACS numbers: 03.75.Ss, 05.30.Fk, 32.80.Pj, 34.50.-s 

By applying a magnetic field to a gas of ultra-cold 
atoms, it is possible to tune the strength and the sign 
of the effective interaction between particles. This phe- 
nomenon, known as Feshbach resonance, offers in the case 
of fermions the unique possibility to study the crossover 
between situations governed by Bose-Einstein and Fermi- 
Dirac statistics. Indeed, when the scattering length a 
characterizing the 2-body interaction at low tempera- 
ture is positive, the atoms are known to pair in a bound 
molecular state. When the temperature is low enough, 
these bosonic dimers can form a Bose-Einstein conden- 
sate (BEC) as observed very recently in 40 K and 6 Li 
0, Q . On the side of the resonance where a is negative, 
one expects the well known Bardeen-Cooper-Schrieffer 
(BCS) model for superconductivity to be valid. However, 
this simple picture of a BEC phase on one side of the res- 
onance and a BCS phase on the other is valid only for 
small atom density n. When n\a\ 3 > 1 the system enters 
a strongly interacting regime that represents a challenge 
for many-body theories [4], |5j,|6j an d that now begins to 
be accessible to experiments 

In this letter, we report on Bose-Einstein condensa- 
tion of 6 Li dimers in a crossed optical dipole trap and 
a study of the BEC-BCS crossover region. Unlike all 
previous observations of molecular BEC made in single 
beam dipole traps with very elongated geometries, our 
condensates are formed in nearly isotropic traps. Analyz- 
ing free expansions of pure condensates with up to 4 x 10 4 
molecules, we measure the molecule-molecule scattering 
length a m = 170tgo° nm at a magnetic field of 770 gauss. 
This measurement is in good agreement with the value 
deduced from the resonance position || and the relation 
a m = 0.6 a of ref. [10], Combined with tight confinement, 
these large scattering lengths lead to a regime of strong 
interactions where the chemical potential fi is on the or- 
der of ksTc where Tc — 1.5 fiK is the condensation tem- 
perature. As a consequence, we find an important mod- 
ification of the thermal cloud time of flight expansion 
induced by the large condensate mean field. Moreover, 
the gas parameter n m a^ is no longer small but on the 



order of 0.3. In this regime, the validity of mean field 
theory becomes questionable 0, UJ. We show, m 
particular, that the anisotropy and gas energy released 



during expansion varies monotonically across the Fesh- 
bach resonance. 

Our experimental setup has been described previously 
[Til fl5j |. A gas of 6 Li atoms is prepared in the abso- 
lute ground state 1 1/2, 1/2) in a Nd-YAG crossed beam 
optical dipole trap. The horizontal beam (resp. vertical) 
propagates along x (y), has a maximum power of = 
2W (P% = 3.3 W) and a waist of - 25 fan (~ 40 /mi). 
At full power, the 6 Li trap oscillation frequencies are 
lu x /2it = 2.4(2) kHz, u y /2n = 5.0(3) kHz, and uj z /2tt = 
5.5(4) kHz, as measured by parametric excitation, and 
the trap depth is ~ 80 fiK. After sweeping the mag- 
netic field B from 5 G to 1060 G, we drive the Zee- 
man transition between 1 1/2, 1/ 2) and 1 1/2, — 1/2) with 
a 76 MHz RF field to prepare a balanced mixture of the 
two states. As measured very recently the Fesh- 
bach resonance between these two states is peaked at 
822(3) G, and for B=1060G, a = -167nm. After 100 ms 
the coherence between the two states is lost and plain 
evaporation provides Nf = N± = N tot /2 = 1.5 x 10 5 
atoms at 10/iK=0.8T F , where k B T F = h 2 kl/2m = 
) 1/3 = ?iw(3iVt t) 1/3 is the Fermi energy. 
Lowering the intensity of the trapping laser to 0.1 Po, the 
Fermi gas is evaporatively cooled to temperatures T at 
or below 0.2 T-p and AT tot w 7 x 10 4 . 

Then, sweeping the magnetic field to 770 G in 200 ms, 
the Feshbach resonance is slowly crossed. In this pro- 
cess atoms are reversibly transformed into cold molecules 
[14L [lfj near the BEC critical temperature as presented 
in figure la. The onset of condensation is revealed by bi- 
modal and anisotropic momentum distributions in time 
of flight expansions of the molecular gas. These images 
are recorded as follows. At a fixed magnetic field, the 
optical trap is first switched off. The cloud expands typ- 
ically for 1 ms and then the magnetic field is increased 
by 100 G in 50 /is. This converts the molecules back into 
free atoms above resonance without releasing their bind- 
ing energy |3J. Switching the field abruptly off in 10 /is, 
we detect free 6 Li atoms by light absorption near the D2 
line. Using this method, expansion images are not al- 
tered by the adiabatic following of the molecular state to 
a deeper bound state during switch-off as observed in our 
previous work |l4j |. Furthermore, we check that there are 
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no unpaired atoms before expansion. In figure lb, a Bose- 
Einstein condensate of 7 Li atoms produced in the same 
optical trap is presented. The comparison between the 
condensate sizes after expansion reveals that the mean 
field interaction and scattering length are much larger 
for 6 L\2 dimers (Fig. la) than for 7 Li atoms (Fig. lb). 

To measure the molecule-molecule scattering length, 
we produce pure molecular condensates by taking ad- 
vantage of our crossed dipole trap. We recompress the 
horizontal beam to full power while keeping the verti- 
cal beam at the low power of 0.035 Pq corresponding 
to a trap depth for molecules U = 5.6 fiK. Tempera- 
ture is then limited to T < 0.9 fiK assuming a conser- 
vative 7] = U /ksT = 6, whereas the critical tempera- 
ture increases with the mean oscillation frequency. Con- 
sequently, with an axial (resp. radial) trap frequency 
of 440 Hz (resp. 5 kHz), we obtain T/T° < 0.3, where 
Tg = huj(0.82N tot /2) 1 / 3 =2.7 /iK is the non interacting 
BEC critical temperature. Thus, the condensate should 
be pure as confirmed by our images. After 1.2 ms of 
expansion, the radius of the condensate in the x (resp. 
y) direction is R x = 51 fim (R y = 103 fim). The re- 
sulting anisotropy R y /R x — 2.0(1) is consistent with the 
value 1.98 0] predicted the scaling equations [lil llfSj. 
Moreover, this set of equation leads to an in-trap radius 
R x = 26/xm (resp. R y — 2.75/zm). We then deduce the 
molecule-molecule scattering length from the Thomas- 
Fermi formula R xy = a ho Q/uj x ^(l5N tot a IIL /2ai lo ) 1 ^, 

with dho — y/h/2mui. Averaging over several images, 
this yields a m =170i 60 nm at 770 G,. Here, the statisti- 
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FIG. 1: Onset of Bose- Einstein condensation in a cloud of 2 x 
10 4 6 Li dimers at 770 G (a) and of 2 x 10 4 7 Li atoms at 610 G 
(b) in the same optical trap, (a): 1.2ms expansion profiles 
along the weak direction x of confinement, (b): 1.4 ms expan- 
sion. The different sizes of the condensates reflect the large 
difference in scattering length a m = 170 nm for 6 Li dimers and 
at — 0.55 nm for 7 Li. Solid line: Gaussian+Thomas-Fermi 
fit. Dashed line: gaussian component. Condensate fractions 
are both 28%. lo x /2tv = 0.59(4) kHz, u y /2ir = 1.6(1) kHz, 
and lu z /2ty = 1.7(1) kHz in (a). lo x /2tt = 0.55(4) kHz, 
lo v /2tt = 1.5(1) kHz, and lj z /2tv = 1.6(1) kHz in (b). 



cal uncertainty is negligible compared to the systematic 
uncertainty due to the calibration of our atom number. 
At this field, we calculate an atomic scattering length of 
a = 306 nm. Combined with the prediction a m = 0.6 a of 
(lpf. we obtain a m = 183 nm in good agreement with our 
measurement. For 7 Li, we obtain with the same analysis 
a much smaller scattering length of O7=0. 65(10) nm at 
610 G also in agreement with theory 20]. 

Such large values of a m bring our molecular conden- 
sates into a novel regime where the gas parameter Jimd^ 
is no longer very small. Indeed, a m = 170 nm and 
n m = 6 x 10 13 cm~ 3 yield n m af n = 0.3. As a first 
conseque nce, corrections due to beyond mean field ef- 
fects |ll|, |21| or to the underlying fermionic nature of 
atoms may play a role, since the average spacing between 
molecules is then of the order of the molecule size ~ a/2. 
Second, even in a mean field approach, thermodynamics 
is expected to be modified. For instance, in the condi- 
tions of Fig. la, we expect a large shift of the BEC criti- 
cal temperature [ill Iha fl3| . The shift calculated to first 
order in rt 1/3 a ATq/TQ = —1.4, is clearl y in applica- 
ble and a more refined approach is required |22|. Third, 
we observe that partially condensed cloud expansions are 
modified by interactions. Indeed, double structure fits 
lead to temperatures inconsistent with the presence of a 
condensate. In Fig. 1, we find T = 1.6 fiK, to be com- 
pared with Tq = 1.4 /iK, whereas for the 7 Li condensate 
T = 0.7/iK=0.6T,S. 

This inconsistency results from the large mean field 
interaction which modifies the thermal cloud expansion. 
To get a better estimate of the temperature, we rely 
on a release energy calculation. We calculate the Bose 
distribution of thermal atoms in a mexican hat poten- 
tial that is the sum of the external potential and the 
repulsive mean field potential created by the conden- 
sate. For simplicity we neglect the mean field result- 
ing from the thermal component. The release energy is 
the sum of the thermal kinetic energy, condensate in- 
teraction energy, and Hartree-Fock interaction energy 
between the condensate and thermal cloud. The tem- 
perature and chemical potential are then adjusted to 
fit the measured atom number and release energy. For 
figure la, we obtain a condensate fraction of 28 % and 
(j, = hui/2(15N tot a m /2a ho ) 2 / 5 = 0.45 //K. The tempera- 
ture T = 1.1 [iK is then found below T° = 1.4 ^K. 

The condensate lifetime is typically ~300ms at 715 G 
(a m = 66 nm) and ~3s at 770 G (a m = 170 nm), whereas 
for a = —167 nm at 1060 G, the lifetime exceeds 30 s. 
On the BEC side, the molecule-molecule loss rate con- 
stant is G = 0.26to;o6 x 10" 13 cm 3 /s at 770 G and 
G = 1.75±g;| x 10" 13 cm 3 /s at 715 G with the fit proce- 
dure for condensates described in [2^] . Combining similar 
results for four values of the magnetic field ranging from 
700 G to 770 G, we find G cx a- 19±0 - 8 . Our data are in 
agreement with the theoretical prediction G cx a~ 2 - 55 of 
ref. [Toj ] and with previous measurements of G in a ther- 
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mal gas at 690 G [ij or in a BEC at 764 G 8]. A similar 
power law was also found for 40 K [iij] . 

We now present an investigation of the crossover from 
a Bosc-Einstcin condensate to an interacting Fermi gas 
(Fig. 2 and 3). We prepare a nearly pure condensate with 
3.5 xlO 4 molecules at 770 G and recompress the trap to 
frequencies of ui x = 2tt x 830 Hz, uj v = 2ir x 2.4 kHz, and 
uj z = 2ir x 2.5 kHz. The magnetic field is then slowly 
swept at a rate of 2 G /ms to various values across the 
Feshbach resonance. The 2D momentum distribution af- 
ter a time of flight expansion of 1.4 ms is then detected 
as previously. 
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FIG. 2: Integrated density profiles across the BEC-BCS 
crossover region. 1.4 ms time of flight expansion in the ax- 
ial (a) and radial (b) direction. The magnetic field is varied 
over the whole region of the Feshbach resonance from a > 
(B < 810 G) to a < (B > 810 G). •: Feshbach resonance 
peak. 



Fig. 2 presents the observed profiles (integrated over 
the orthogonal direction) for different values of the mag- 
netic field. At the lowest field values B < 750 G, 
n>m.a>m ^ 1) condensates number are relatively low be- 
cause of the limited molecule lifetime. As B increases, the 
condensate width gradually increases towards the width 
of a non interacting Fermi gas, and nothing dramatic 
happens on resonance. At the highest fields (B> 925 G), 
where k-p\a\ < 3, distributions are best fitted with zero 
temperature Fermi profiles. More quantitatively, Fig. 3b 
presents both the gas energy released after expansion E rc \ 
and the anisotropy 77 across resonance. These are calcu- 
lated from gaussian fits to the density after time of flight: 
E re i = m(2(jy + a x )/2r 2 and 77 = a y /a x , where Oi is the 
rms width along i, and r is the time of flight ^t|- O n the 
BEC side at 730 G, the measured anisotropy is 77 ~1.6(1) 
, in agreement with the hydrodynamic prediction, 1.75. 
It then decreases monotonically to 1.1 at 1060 G on the 
BCS side. On resonance, at zero temperature, superfluid 
hydrodynamic expansion is expected [25| corresponding 



to 77 — 1.7. We find however 77 = 1.35(5), indicating a 
partially hydrodynamic behavior that could be due to a 
reduced superfluid fraction. On the a < side, the de- 
creasing anisotropy would indicate a further decrease of 
the superfluid fraction that could correspond to the re- 
duction of the condensed fraction of fermionic atom pairs 
away from resonance observed m_0> @ • Interestingly, our 
results differ from that of ref.|26J where hydrodynamic 
expansion was observed at 910 G in a more elongated 
trap for T/T F ~ 0.1. 
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FIG. 3: (a): scattering length between the 1 1/2, 1/2) and 
|l/2, —1/2) 6 Li states. The Feshbach resonance peak is lo- 
cated at 820 G (dotted line), (b): anisotropy of the cloud, (c): 
release energy across the BEC-BCS crossover region. In (c), 
the dot-dashed line corresponds to a T = ideal Fermi gas. 
The dashed curve is the release energy from a pure conden- 
sate in the Thomas- Fermi limit. The solid curve corresponds 
to a finite temperature mean field model described in the text 
with T = 0.5 Tg. Arrow: fc P |a| = 3. 



In the BEC-BCS crossover regime, the gas energy re- 
leased after expansion E rc \ is also smooth (Fig. 3c). E IC \ 
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presents a plateau for B < 750 G, and then increases 
monotonically towards that of a weakly interacting Fermi 
gas. The plateau is not reproduced by the mean field 
approach of a pure condensate (dashed line). This is a 
signature that the gas is not at T = 0. It can be under- 
stood with the mean field approach we used previously 
to describe the behavior of the thermal cloud. Since the 
magnetic field sweep is slow compared to the gas col- 
lision rate we assume that this sweep is adiabatic 
and conserves entropy [U We then adjust this entropy 
to reproduce the release energy at a particular magnetic 
field, B — 720 G. The resulting curve as a function of 
B (solid line in Fig. 3c) agrees well with our data in the 
range 680 G< B < 800 G, where the condensate frac- 
tion is 70%, and the temperature is T w T<S/2 = 1.2 fiK. 
This model is limited to n m a^, < 1. Near resonance 
the calculated release energy diverges and clearly departs 
from the data. On the BCS side, the release energy of 
a T = ideal Fermi gas gives an upper bound for the 
data (dot-dashed curve), as expected from negative in- 
teraction energy and a very cold sample. This low tem- 
perature is supported by our measurements on the BEC 
side and the assumption of entropy conservation through 
resonance which predicts T = 0.1 Tp |27j . 

On resonance the gas is expected to reach a univer- 
sal behavior, as the scattering length a is not a rele- 
vant parameter any more In this regime, the re- 
lease energy scales as E Te \ — \/l + ]3~E® el , where E° el is 
the release energy of the non-interacting gas and (3 is a 
universal parameter. From our data at 820 G, we get 
f3 = —0.64(15) . This value is larger than the Duke re- 
sult ft = -0.26 ±0.07 at 910 G [Iff, but agrees with that 
of Innsbruck (3 = -0.68±£;^ at 850 G @, and with the 
most recent theoretical prediction = —0.56 0. Around 
925G, where a = -270 nm and (fcpH)- 1 = 0.35, the re- 
lease energy curve displays a change of slope. This is 
a signature of the transition between the strongly and 
weakly interacting regimes. It is also observed near the 
same field in Q through in situ measurement of the 
trapped cloud size. Interestingly, the onset of resonance 
condensation of fermionic atom pairs observed in 40 K Q 
and 6 Li [{J, corresponds to a similar value of fcp|a|. 

In summary, we have explored the whole region of the 
6 Li Feshbach resonance, from a Bose-Einstein condensate 
of fermion dimers to an ultra-cold interacting Fermi gas. 
The extremely large scattering length between molecules, 
that we have measured leads to novel BEC conditions. 
We have observed hydrodynamic expansions on the BEC 
side and non-hydrodynamic expansions at and above res- 
onance. We suggest that this effect results from a reduc- 
tion of the superfluid fraction and we point to the need of 
a better understanding of the dynamics of an expanding 
Fermi gas. 
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